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CALCULATION  OF  FAR  FIELD,  SHIP  RADIATION 
AND  DIFFRACTION  WAVES 


1.  INTRODUCTION 

This  memorandum  report  outlines  the  initial  set  of  proposed  approaches  to  obtain  the  far  field 
radiation  and  diffraction  waves  generated  by  a  ship  advancing  at  constant  speed  and  oscillating  in  the 
presence  of  ambient  ocean  waves.  Radiation  waves  are  caused  by  the  oscillation  of  the  ship  while  the 
diffraction  waves  refer  to  the  disturbance  to  the  incoming  ambient  waves  caused  by  the  presence  of  the 
ship. 


The  report  first  presents  the  linearized  mathematical  formulation  which  is  often  used  to  solve 
these  problems.  The  assumptions  underlying  this  formulation  are  briefly  stated.  The  further 
simplifications  which  are  required  to  arrive  at  the  commonly  used  strip  theory  approach,  as  imple¬ 
mented  in  the  DTNSRDC  Ship  Motion  Program,  are  described.  It  is  pointed  out  that  this  is  a  two- 
dimensional  near  field  solution  which  is  accurate  near  the  ship  hull.  The  mathematical  convenience  as 
well  as  the  limitations  resulting  from  these  simplifications  are  discussed.  The  principal  differences  and 
similarities  of  the  radiation  and  diffraction  problems  in  the  strip  theory  approach  are  indicated. 

The  principal  interest  in  the  large  majority  of  previous  studies  has  been  the  calculation  of  the  near 
field  flow  in  order  to  obtain  the  added  mass,  damping,  and  exciting  forces  on  the  ship.  The  two  steps 
required  to  convert  the  near  field  solutions  to  the  calculation  of  the  far  field  wave  pattern  are  outlined 
in  some  detail.  The  first  step  is  to  convert  the  near  field  solution  to  an  equivalent  distribution  of  far 
field  singularities  placed  on  or  near  the  hull  of  the  ship.  The  second  step  requires  the  asymptotic 
evaluation  of  the  far  field  wave  pattern  due  to  these  singularities. 

Three  approaches  are  described  for  obtaining  the  far  field  singularity  distribution.  The  first 
approach  involves  directly  converting  the  calculated  near  field  two-dimensional  singularity  distribution 
to  an  equivalent  distribution  of  three-dimensional  far  field  singularities  by  using  the  method  of  matched 
asymptotic  expansions.  The  second  approach  is  somewhat  similar  to  the  first.  It  involves  first  convert¬ 
ing  the  calculated  or  measured  forces  on  the  two-dimensional  strips  to  a  distribution  of  near  field  singu¬ 
larities  and  then  using  the  first  approach  to  obtain  the  far  field  singularity  distribution.  This  approach 
makes  use  of  the  abundant  data  on  the  added  mass,  damping,  and  exciting  forces  acting  on  ship  hulls. 
The  third  approach  characterizes  the  entire  singularity  distribution  on  the  ship  hull  by  an  amplitude 
function  which  represents  an  integrated  effect  of  the  distribution. 

It  is  pointed  out  that  asymptotic  expressions  have  been  derived  for  the  far  field  behavior  of  pulsat¬ 
ing  sources  and  dipoles  which  are  otherwise  stationary.  Thus,  the  evaluation  of  the  far  field  wave  pat¬ 
tern  for  ships  with  zero  speed  is  relatively  simple.  The  situation  is  quite  different  for  the  case  of  for¬ 
ward  speed,  for  which  the  formulas  and  wave  pattern  of  the  singularities  are  considerably  more  com¬ 
plex.  There  have  been  relatively  few  studies  which  consider  the  far  field  behavior  of  these  singularities. 
These  studies  are  concerned  with  obtaining  some  general  features  of  the  far  field  behavior  and  are 
largely  written  in  French  and  German.  In  view  of  the  complexity  of  the  forward  speed  case,  some 
alternate  simpler  approaches  (which  may  be  suitable  for  low  and  moderate  speed  cases)  are  briefly  indi¬ 
cated. 
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2.  LINEARIZED  FORMULATION  OF  PROBLEM 


The  problem  is  most  often  formulated  for  a  coordinate  system  (x,  y,  z),  shown  in  Fig.  1,  moving 
in  the  x-direction  with  the  mean  forward  speed  U  of  the  ship.  The  moving  coordinate  system  (x,  y,  z) 
is  related  to  the  fixed  spatial  coordinate  system  (x0,  yo<  zo)  by 


(x,  y,  z)  =  (x0-  Ut,  y0,  z0)  (1) 

where  t  is  time.  Also,  the  frequency  of  encounter  w  in  the  moving  coordinate  system  is  related  to  the 
frequency  a>0  of  incident  waves  in  the  fixed  coordinate  system  by 


a >  =  o>0  —  Uko  cos  (3  =  co0  — 


Uajj 


8 


where  /3  is  the  angle  between  the  direction  of  propagation  of  the  incident  wave  and  the  x-direction 
ko  =  wo/ 8  =  27t/\  is  the  wavenumber 

A.  is  the  wavelength. 


(2) 


The  assumptions  are  made  that  the  fluid  is  inviscid  and  incompressible,  the  flow  is  irrotational, 
and  surface  tension  effects  may  be  neglected.  Under  these  assumptions,  the  flow  field  is  given  by  the 
velocity  potential  (x,  y,  z)  which  satisfies  Laplace’s  equation 


V2<6(x,  y,  z) 


82<P  92<t> 

8x2  by2 


(3) 


The  total  potential  <I>  may  be  written  as  the  sum  of  the  steady  potential  £70,  which  gives  rise  to  the 
Kelvin  waves,  and  the  unsteady  potential  6 


T>  =  U<f>  +  (4) 

The  principal  interest  in  the  present  report  is  the  unsteady  potential  <j>.  Assuming  that  the  amplitudes 
of  the  incident,  radiation,  diffraction,  and  Kelvin  waves  are  all  small  leads  to  the  linearized  formulation 
given  below.  Nonlinear  formulations  are  given  by  Ogilvie  and  Tuck,1  Newman,2  and  Maruo.3  How¬ 
ever,  it  should  be  noted  that  even  the  linearized  problem  is  of  considerable  difficulty,  and  most  of  the 
current  investigations  focus  on  solving  various  simplifications  of  the  linearized  formulation. 


Under  the  above  linearizing  assumptions,  the  potential  6  may  be  further  decomposed  into 
separate  potentials  due  to  the  incident  wave,  each  of  the  six  rigid  body  oscillations  shown  in  Fig.  1,  and 
the  diffracted  wave.  Assuming  the  motions  to  be  sinusoidal  at  the  encounter  frequency  to,  the  notation 
shown  in  Fig.  1  suggests  that  the  motions  of  the  ship  are  given  by 

(fi.  *2.  *3.  fli.  n2,  a 3) <?'"'=  (I,,  £2,  £3,  £4,  fj,  £6')ey“'.  (5) 

The  unsteady  potential  4>  may  then  be  written  as 

4>  =  [A  (<£0  +  fa)  +  ]T  t‘j<t>j]e,al  (6a) 

j- 1 

where  A  is  the  amplitude  of  the  incident  wave, 

<f> 0  =  —  exp  [/c0(z  -  ix  cos  /3  —  iy  sin  /!)],  (6b) 

G>0 

is  the  potential  for  an  incident  wave  of  unit  amplitude  whose  direction  of  propagation  is  at  angle  /3  with 
the  x-axis,  and  7  is  the  diffraction  potential.  The  original  problem  is  thus  restated  as  the  separate  solu¬ 
tion  of  the  diffraction  problem  where  incident  waves  act  upon  the  equilibrium  submerged  hull  surface, 
and  the  radiation  problem  where  the  ship  undergoes  prescribed  oscillations  in  the  absence  of  incident 
waves. 
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On  the  mean  free  surface  2  =  0  each  of  the  above  potentials  0  ,  j  =  0  to  7,  must  satisfy  the 
linearized  free  surface  condition 


dt2 


320 


Ut,  y,  z)  +  g-— = 


-  U— 
dt  dx 


4>j  +  g 


dz 

3z 


2U 


30 ,  30 


9r2  ~~  3r  9.v 

on  z  =  0,  j  =  0,  1,  . . . ,  6,  7. 


+  (72 


i  9x6  ,• 
3  .v' 


+  s- 


30; 

9z 


=  0 


(7) 


On  the  mean  submerged  hull  surface  5,  the  diffraction  potential  07  must  satisfy  the  kinematic  condition 


307 

9/2 


300 

9/2 


on  5 


(8) 


where  n  is  the  normal  to  the  hull  surface.  The  kinematic  condition  on  S  to  be  satisfied  by  the  radiation 
potentials  for  the  case  of  a  ship  moving  at  forward  speed  U  have  been  derived  by  Timman  and  New¬ 
man4 


tL 1 

9/2 


=  {ioja  +  Vx  [a  x  V(6r0)]}  •  n 


(9) 


where  a  is  the  displacement  of  a  point  on  S.  Recalling  that  the  disturbance  associated  with  the  steady 
forward  motion  is  small,  and  expressing  a  in  terms  of  the  oscillatory  translations  J  and  rotations  O 
defined  in  Eq.  (5)  and  Fig.  1,  the  above  equation  becomes2:5 


3  0; 

dn 


=  icorij  +  Urrij  on  S,  j  =  1,  . , . ,  6 


(10) 


where 

> 

(n\,  /22‘,  223)  =  77  (11) 

(/24,  n5,  n6)  =  (7  x  77)  (12) 


/■  is  the  radius  vector  from  the  origin  to  a  point  on  the  hull  surface 
(m i,  m2,  nis,  2224)  =>0 

//25  =  n3,  m6  =  -  «2.  (13) 

The  solution  of  even  the  above  linearized  problem  is  a  formidable  task.  For  example,  the  com¬ 
puter  cost  for  the  method  developed  by  Chang6  to  solve  the  above  problem  is  nearly  two  orders  of 
magnitude  larger’  than  the  corresponding  cost  for  the  widely  used  DTNSRDC  Ship  Motion  Program 
(SMP)3'8,  which  uses  the  simpler  strip  theory  approach.  This  theory,  which  is  described  in  greater 
detail  later  in  this  report,  essentially  converts  the  complex  three-dimensional  problem  to  a  series  of 
simpler  two-dimensional  problems. 

3.  SIMPLIFICATIONS  TO  LINEARIZED  THEORY 

Most  of  the  simplifications  to  the  above  linearized  theory  are  based  on  making  one  or  more  of  the 
following  assumptions  regarding  the  magnitude  of  the  ship  and  wave,  parameters2, 3,9 

B/L  «  1  (14a) 
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T/L  «  1 

(14b) 

co2L/g  =  2-nL/k  «  1 

(14c) 

uL/g  =  2ttL/K  »  1 

(14d) 

where  B  is  the  beam  (width)  of  the  ship 

L  is  the  length  of  the  ship 

Tis  the  draft  of  the  ship. 

Assumptions  (14a)  and  (14b)  lead  respectively  to  thin  ship  and  flat  ship  (or  planing  ship)  theories. 
Assumptions  (14c)  and  (14d)  lead  respectively  to  low  frequency  (or  large  wavelength)  and  high  fre- 
quency  (or  small  wavelength)  theories.  A  relatively  large  number  of  studies  are  based  on  slender  body 
theory,  in  which  assumptions  (14a)  and  (14b)  are  concurrently  used.  A  particular  form  of  slender  body 
theory,  in  which  assumption  (14d)  is  also  used,  leads  to  the  widely  used  strip  theory  which  is  described 
in  detail  in  the  following  section. 


For  motions  in  the  vertical  plane  (surge  fj,  heave  £3,  and  pitch  £5),  thin  ship  theory  leads  to  a 
simple,  direct  solution  for  the  potential  0,  as  follows.  By  using  Green’s  Theorem,  0  can  generally  be 
expressed  as  the  following  integral  over  the  hull  surface  S 


*-tss 


r  ^0 
(js  —  0 

dn 


SGc 


dn 


c/S 


(15) 


where  Gs  is  a  pulsating  source  which  satisfies  Laplace’s  Eq.  (3)  and  the  free  surface  condition  (7). 
The  use  of  thin  ship  theory  allows  the  integral  to  be  evaluated  over  the  vertical  (y  =  0)  centerplane  of 
the  ship.  The  derivative  d/d n  may  then  be  approximated  by  d/dy.  The  crucial  simplification  results 
from  the  fact  the  that  dGs/dn  =  dGs/dy  =  0  on  the  centerplane  y  =  0  since  Gs  is  an  even  function  in 
v.  Equation  (15)  then  simplifies  to 


W*  (t6> 

where  d<j>/dy  is  the  prescribed  velocity  on  the  ship  hull  and  is  given  in  Eqs.  (8)  to  (13).  The  above 
equation  then  represents  a  direct  solution  for  0. 


For  the  case  ot  a  nonoscillating  ship  advancing  at  constant  speed  U,  thin  ship  theory  yields  useful 
results  for  the  steady-state  flow  field.  For  the  case  of  a  ship  executing  oscillations  in  the  vertical  plane, 
Peters  and  Stoker  obtain  the  disappointing  result  that  the  oscillations  cause  no  disturbance  to  the 
fluid,  i.e.,  the  added  mass  and  damping  forces  are  zero.  Mathematically,  the  difference  between  the 
steady-state  and  oscillating  cases  may  be  viewed  as  follows.  The  thin  ship  assumption,  given  by  Eq. 
(14a),  may  be  rewritten  as 

B/L  =  0  (es),  eB  «  1.  (17) 

For  a  ship  advancing  at  finite  velocity,  the  following  dimensionless  form  for  U  may  be  taken  to  be 
order  1 


F2=  U2/gL  =  0(1)  (18) 

where  F  is  the  Froude  number.  The  linearizing  assumptions  stated  under  Eq.  (4)  imply  that  the  ship 
oscillations  t/j  are  small 

f/  =  0  («*>.  «  1.7=1,  • . 6.  (19) 

Thus,  the  flow  disturbance  caused  by  the  steady  velocity  is  of  first  order,  while  the  disturbance  caused 
by  the  ship  oscillations  is  of  higher  order.  Physically,  the  results  of  Peters  and  Stoker  simply  state  that 
a  thin  vertical  disk  excuting  small  oscillations  in  the  vertical  plane  creates  disturbances  which  are  negli¬ 
gible  to  first  order.  Newman  has  carried  out  the  analysis  to  second  order  and  obtained  nonzero 
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expressions  for  the  added  mass  and  damping  forces.  However,  to  quote  Newman,  "Unfortunately,  the 
second-order  equations  are  rather  complex.. .as  might  be  expected  from  the  use  of  a  systematic  pertur¬ 
bation  procedure."  Accordingly,  his  proposed  approach  has  not  gained  acceptance  as  a  practical  compu¬ 
tational  method. 


It  should  be  noted  that  thin  ship  theory  predicts  nonzero  hydrodynamic  disturbances  for  ship 
oscillations  in  the  lateral  y-direction  (sway  |2,  roll  £4,  and  yaw  f6).  However,  the  crucial  direct  solution 
feature  for  motions  in  the  vertical  plane  is  lost  in  this  case.  For  motions  in  the  lateral  direction,  the 
fundamental  singularity  is  the  dipole  GL  given  by 


Gl  = 


dGs 

9y 


(20) 


The  quantity  dGL/dn  =  bGjdy  =  d2Gs/dy2  is  an  even  function  in  y  and  is  not  necessarily  equal  to 
zero  on  the  vertical  centerplane  y  =  0.  Thus  both  terms  on  the  right  hand  side  of  Eq.  (15)  must  be 
retained. 


The  lack  of  success  of  first  order  thin  ship  theory  led  Peters  and  Stoker10  to  consider  flat  ship 
theory,  which  uses  assumption  (14b).  Newman11,12  points  that  this  theory  leads  to  the  considerable 
complexity  of  solving  singular  integral  equations. 


The  low  frequency  assumption  (14c)  leads  to  two  cases  of  general  interest  depending  on  the  mag¬ 
nitude  of  the  Froude  number  F  defined  in  Eq.  (18): 


F2  =  0(1) 
F2  —  0. 


Assumptions  (14c)  and  (21a)  reduce  the  free  surface  condition  given  in  Eq.  (7)  to 


U 2 


9z 


0  on  r  = 


0 


(21a) 

(21b) 

(22) 


and  the  hull  boundary  condition  (10)  to 


90/ 

-r—  =*  U  m,  on  S,  j  *=  1,  . . . ,  6. 
an  J 


(23) 


Thus,  all  dependence  on  a>  is  removed,  resulting  in  a  purely  steady-state  problem  for  a  nonoscillating 
ship.  Assumptions  (14c)  and  (21b)  reduce  the  free  surface  condition  to 


and  the  hull  boundary  condition  to 


(24) 


90,- 

— —  =  itonj  on  S,  j  =  1,  . . .  ,  6.  (25) 

an 


Equation  (25)  shows  that  the  oscillatory  character  of  the  ship  motions  has  been  preserved.  This  large 
wavelength,  low  speed  theory  leads  to  considerable  simplifications.  For  example,  Newman13  shows  that 
the  diffraction  potential  <f>-,  can  be  approximated  by 


90p 

dx 


<t>\  + 


9</>o 

9y 


(f)  2 


90° 

"aT  03 


(26) 


One  simplification  in  the  derivation  of  the  above  equation  is  that  the  derivatives  of  0O  can  be  taken  to 
be  constants,  evaluated  at  a  convenient  location  around  the  body.  The  present  author14  has  exploited 
the  fact  that  this  theory  leads  to  body  translational  motions  which  tend  to  follow  those  of  the  surface 
wave  to  derive  expressions  for  the  forces  acting  on  wide  classes  of  buoys.  However,  this  case  is  not  of 
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primary  interest  for  the  present  study  due  to  the  fact  that  Eq.  (24)  shows  that  the  free  surface  behaves 
like  a  flat  plane.  This  implies  the  absence  of  radiation  waves  (and  the  damping  forces  which  give  rise  to 
these  waves). 

4.  STRIP  THEORY  (NEAR  FIELD  SOLUTION) 

Strip  theory  is  a  particular  form  of  slender  body  theory  which  seeks  to  reduce  the  three- 
dimensional  problem  defined  by  Eqs.  (1)  to  (13)  to  a  series  of  simpler  two-dimensional  problems  for 
cross  sections  (or  strips)  at  various  stations  along  the  longitudinal  x-axis.  The  resultant  theory  gives  a 
near  field  solution  applicable  only  in  the  neighborhood  of  the  hull.  The  extension  to  the  far  field  will 
be  discussed  in  the  following  section.  The  near  field  solution  is  appropriate  for  computing  the  wave 
pattern  near  the  ship  and  the  forces  acting  on  the  hull.  A  number  of  different  versions  of  strip  theory 
and  slender  body  theory  have  been  used  to  analyze  radiation  and/or  diffraction  problems;  see,  for 
example,  the  survey  by  Newman.9  The  description  given  below  of  the  strip  theory  for  the'  radiation 
problem  follows  closely  the  formulation  in  the  widely  used  DTNSRDC  Ship  Motion  Program  (SMP).5,8 
The  description  of  the  solution  for  the  diffraction  potential  differs  from  the  formulation  given  in  SMP. 
In  SMP,  the  diffraction  potential  is  not  computed  explicitly.  Instead,  the  diffraction  forces  are  com¬ 
puted  from  the  radiation  potentials. 

4.1  Radiation  Problem 


In  the  SMP  formulation  for  the  radiation  problem,  the  slender  body  assumptions,  Eqs.  (14a)  and 
(14b),  are  rewritten  in  the  following  form  involving  the  derivatives  and  normals  in  the  neighborhood  of 
the  hull 


4- -0(1),  =  0(l/e), 

ox  dy 


dz 


=  0(1/ «) 


nx  ~  0(1),  ny  =  0(1/ e),  n.  =  0 ( 1/e ) .  (27) 

These  assumptions  reduce  the  three-dimensional  Laplace’s  Eq.  (3)  to  the  following  two-dimensional 
form  in  the  y-z  plane 


32d>; 

3/ 


3 -<t>j 
dz2 


=  0,  j  =  1,  ....  6. 


(28) 


In  order  to  reduce  the  free  surface  condition,  Eq.  (7),  to  two-dimensional  form  near  the  hull,  assump¬ 
tion  (14c)  is  used  in  the  form 


_3_ 

3/ 


»  U  f-. 

ox 


This  reduces  Eq.  (7)  to  the  following  standard  speed  independent  two-dimensional  form 

Hj 


-a>20,  +  g- 


dz 


=  0,  j-  1, 


6. 


(29) 


(30) 


The  speed  dependent  boundary  condition  in  Eq.  (10)  is  eliminated  as  follows.  The  potentials  <f>  are 
divided  into  two  parts. 


4>j-<t>j  +  J=  1 . 6  (31) 

lOJ 

where  the  superscripts  0  and  U  respectively  denote  speed  independent  and  speed  dependent  parts, 
resulting  in  the  revised  boundary  conditions 

Wj/dn  =  mrtj,  j  =  1,  ...  ,  6  (32a) 

3 <t>ju/dn  -  ioirrij,  j  =  l,  ...  ,  6.  (32b) 
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Use  of  the  relations  given  in  Eq.  (13)  leads  to  the  following  expression  for  </>,■ ,  completely  in  terms  of 
the  speed  independent  potentials  <j>j 

<t>j  =  <t>],  J  =  1,  2,  3,  4 
05  =  05°  +  —  03° 

ICO 

06  =  <t>6  -  y-  <t>2 •  (33) 

1(0 

In  addition,  the  assumptions  of  strip  theory  which  imply  that  there  is  no  interaction  between  cross  sec¬ 
tions  lead  to  the  following  simple  relations  between  the  rotational  motions  pitch  ij  =  5)  and  yaw 
(j  =  6)  and  the  corresponding  translational  motions  heave  ( j  =  3)  and  sway  ( j  =  2) 

05°  =  -  x030,  06°  =  *02°-  (34) 

The  formulation  contained  in  Eqs.  (27)  to  (34)  reduces  the  original  three-dimensional  problem  to 
the  solution  of  the  four  fundamental  two-dimensional  speed  independent  potentials 

0>,  z),j-  1,  2,  3,  4 

at  a  series  of  cross  sections  along  the  hull.  The  Frank  close-fit  source  distribution  method  is  used  to 
solve  this  formulation.  Briefly,  the  method  models  only  the  +y  half  of  a  given  cross  section,  which  is 
assumed  to  be  symmetric  with  respect  to  the  vertical  *  -  z  plane,  as  shown  in  Fig.  2.  The  actual 
continuous  half-cross  section  is  approximated  by  a  series  of  N  straight  line  segments  connecting  N  4-  1 
input  points  (y,-,  z,)  where  y,-  is  the  half-breadth  of  the  section  at  depth  z,-.  A  pulsating  line  source  of 
strength  Q,  is  placed  over  each  line  segment,  where  Q,  varies  from  segment  to  segment.  The  source 

strengths  over  the  —y  half  of  the  section  are  symmetric  for  the  vertical  modes  (j  =•  1,  3)  and  antisym¬ 

metric  for  the  lateral  modes  (j  =  2,  4),  as  follows 

Q-S-y„  4)  =  <2/0/,  a)  7=1.  3  (35a) 

QS-y-n  zi )  =  -Qiiyi,  2,)  j  =  2,  4.  (35b) 

The  kinematic  boundary  conditions  given  in  Eq.  (32a)  are  used  to  obtain  a  system  of  //complex  alge¬ 
braic  equations  for  the  complex  strengths  Qr  This  procedure  is  then  repeated  at  a  series  of  cross  sec¬ 
tions  along  the  length  of  the  ship. 

The  effect  of  forward  speed  is  contained  in  the  expressions  for  05  and  06  given  in  Eq.  (33)  (which 
are  directly  related  to  <f>fi  and  in  the  expression  for  the  radiation  forces 


where  Fj2d  is  the  two-dimensional  radiation  force  acting  on  a  given  cross  section. 

C  is  the  perimeter  of  the  cross  section 

Zk  is  the  amplitude  of  the  Alh  mode  of  motion. 

The  radiation  forces  Fj  on  the  entire  hull  are  then  given  by 

FJ  =  JL  fj2d(x)(Lx,  j  —  l,  ,  6.  (37) 

In  spite  of  the  many  simplifications  listed  above,  SMP  appears  to  give  reasonably  accurate  results 
for  forces  and  motions  compared  to  experiment5  and  to  a  more  complex  three-dimensional  approach.6 
Thus,  it  will  be  used  in  the  present  study  as  the  initial  approach  for  obtaining  the  near  field  solution.  A 
more  complex  unified  strip  theory  developed  by  Newman  and  Sclavounos,16  which  accounts  for  an 
interaction  between  ship  cross  sections,  may  be  used  at  a  later  stage  if  the  additional  accuracy  is 
required. 
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It  should  be  noted  that  there  are  a  number  of  inconsistencies  in  the  above  theory.  Perhaps  the 
most  noticeable  is  the  neglect  of  the  convection  term  Ud/dx'm  the  free  surface  condition,  Eq.  (30) 
and  its  inclusion  in  the  force  Eq.  (36).  Ogilvie  and  Tuck1  have  developed  a  more  consistent  strip 
theory  which  shows  that  most,  but  not  all,  of  the  higher  order  effects  tend  to  cancel  each  other.  The 
effect  of  the  nonzero  higher  order  corrections,  which  affect  only  the  cross-coupling  force  coefficients,  is 
not  yet  well  estblished. 

4.2  Diffraction  Problem 


SMP  does  not  explicitly  calculate  the  diffraction  potential  <f>7.  Instead,  the  diffraction  forces  F° 
are  computed  from  the  values  of  9<£o /dn  (where  d>o  is  given  in  Eq.  (6b))  and  the  previously  calculated 
radiation  potentials  and  <f>/  (Eqs.  (31)-(34))  on  the  hull  surface  by  using  a  generalized  form  of  the 
Haskind  relations,  as  follows 


-Ft" 

ICO 


v J 


dn 


dS,  j=  1,  ...  ,  6. 


(38) 


The  Haskind  relations  were  originally  derived  for  the  case  of  zero  speed.  However,  Newman17  general¬ 
ized  these  relations  to  the  forward  speed  case  for  rather  general  conditions.  The  principal  case  where 
the  relations  do  not  apply  is  for  a  thin  ship  in  head  seas,  /3  =  180  deg. 


The  boundary  condition  18)  for  <f>i  leads  to  several  differences  in  the  formulation  from  the 
previously  discussed  case  for  the  radiation  potentials  <f>j.  Equations  (8)  and  (6b)  show  that 

d<t>i  —dd>o  —9  ig 

=  ~Tn~  =  a7  ^ exp  [k° (z  ~  ix  cos  sin 


”  f(k o  =  <»6 Vs,  X,  y,  z,  /3).  (39) 

Oh  the  other  hand,  the  four  fundamental  radiation  potentials  <£?  j  =  1,  2,  3  4  have  the  functional 
form 


<f>j  =  /(<«>.  y.  z).  (40) 

Since  the  x-dependence  of  <f>i  is  periodic,  it  is  usual  to  write  <bi  as  the  following  product 


4>t(x,  y,  z)  =  z,  ca0,  /3)  e  ,x 


(41) 


where  /0  =  k0  cos  /3 

^7  is  a  two-dimensional  solution. 

Substitution  of  the  slender  body  assumptions  (27),  and  Eq.  (41)  into  the  three-dimensional  Laplace’s 
Eq.  (31  results  in  the  following  two-dimensional  Helmholtz  equation  for  <J>7  {y,  z) 


By2 


(42) 


as  opposed  to  a  two-dimensional  Laplace’s  equation  for  the  radiation  potentials.  A  second  difference  is 
that  the  substitution  of  Eq.  (41)  along  with  assumption  (27)  into  the  free  surface  equation  (7)  yields 


—OJQ  <b7  +  g 


(43) 


whereas  the  free  surface  Eq.  (30)  for  is  in  terms  of  the  encounter  frequency  w. 


The  need  to  solve  the  Helmholtz  Eq.  (42)  instead  of  Laplace’s  Eq.  (28)  makes  the  diffraction 
problem  somewhat  more  difficult  since  the  singularities  appropriate  for  the  Helmholtz  equation  must  be 
used.  Faltinsen,18  Troesch,19  and  Maruo3  solve  the  above  formulation  for  zero  and  low  speed  cases. 
Faltinsen  finds  that  "It  would  be  time  consuming  to  evaluate  the  solutions. ...for  a  ship  with  arbitrary 
cross  sections,  and  obtains  results  for  only  circular  cross  section  cases. 
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Strictly  speaking,  the  dependence  of  07  on  x  shown  in  Eq.  (41)  violates  the  strip  theory  assump¬ 
tion  that  there  is  no  interaction  between  cross  sections.  For  this  reason,  strip  theory  approaches  either 
bypass  the  calculation  of  07,  as  in  the  previously  described  case  of  SMP  which  uses  the  Haskind  rela¬ 
tions  to  calculate  the  diffraction  forces,  or  attempt  to  eliminate  the  x-dependence  shown  in  Eq.  (41)  by 
setting  /0  =  0.  This,  in  effect,  simplifies  the  Helmholtz  Eq.  (42)  to  the  simpler  Laplace’s  Eq.  (28). 
One  obvious  case  where  /0  =  0  is  for  beam  seas,  j8  =  90°.  Another  case  would  appear  to  be  for  long 
waves  such  that 


k0 


2tt 

X 


g 


«  1. 


(44) 


However,  this  would  violate  the  strip  theory  assumption  (14d)  that  the  waves  are  short.  Newman13 
shows  that  strip  theory  is  applicable  to  the  diffraction  problem  over  a  restricted  range  of  X  such  that  it  is 
short  compared  to  the  ship  length  L  but  long  compared  to  the  beam  B 


B  «  X  «  L.  (45) 

In  this  restricted  range,  the  diffraction  potential  is  solved  in  a  manner  which  is  very  similar  to  the  radia¬ 
tion  potentials13 


07 - i(D0  (l//2  — 

302  ,  v 
“3 —  O',  z)  =  ny 

an  * 


i  sin  /303)  elkoXC°50 
303 


0« 


(y,  z)  -  n. 


(46a) 

(46b) 


where  02  and  03  are  respectively  the  antisymmetric  and  symmetric  two-dimensional  potentials  which 
satisfy  Laplace’s  Eq.  (28)  and  the  free  surface  condition  (43).  For  the  case  of  /3  =  45  deg.  Troesch19 
shows  that  the  pressure  distributions  on  the  hull  of  an  ore  carrier  obtained  by  using  this  simpler 
approach  show  only  moderate  differences  from  those  obtained  by  using  the  more  complex  Helmholtz 
Eq.  (42).  Hence,  the  formulation  outlined  in  Eqs.  (46a)  and  (46b)  will  initially  be  used  to  compute  07. 


The  above  shows  that  two-dimensional  solutions  for  the  diffraction  potential  are  more  difficult 
and/or  are  applicable  over  a  more  restricted  range  of  parameters  than  the  corresponding  solutions  for 
the  radiation  potentials.  For  fully  three-dimensional  approaches,  such  as  that  developed  by  Chang,6 
there  is  less  difference  between  the  radiation  and  diffraction  problems  since  both  are  dependent  on  x. 


5.  FAR  FIELD  DISTRIBUTION  OF  SINGULARITIES 


The  near  field  singularity  distribution  obtained  by  the  above  strip  theory  must  be  converted  to  a 
singularity  distribution  which  is  appropriate  for  calculating  the  far  field  wave  pattern.  One  obvious 
shortcoming  of  the  near  field  solution  is  that  the  pulsating  sources  resulting  from  the  two-dimensional 
formulation  gives  waves  of  constant  amplitude  at  large  distances  from  the  hull.  Physically,  the  waves 
attenuate  to  zero  at  large  distances  due  to  cylindrical  spreading. 

The  present  section  will  discuss  three  approaches  for  obtaining  the  equivalent  far  field  singularity 
distribution.  Since  the  formulation  which  will  be  initially  used  for  the  diffraction  potential  07  is  similar 
to  the  formulation  for  the  radiation  potentials  <f>j,  j  =  1,  ...  ,  6,  the  discussion  essentially  applies  to 
both  types  of  potentials. 


Far  from  the  ship,  the  details  of  the  ship  hull  are  lost,  and  the  hull  disturbance  may  be  taken  to 
be  caused  by  a  line  of  pulsating  singularities  placed  on  the  x-axis  along  the  length  L  of  the  ship.  In 
general,  the  outer  potential  07  can  be  expressed  in  the  form 


*/-S, 


Q ’;(£)  +  dj(& 


G3D(x  -  y,  z)  d£ 


(47) 
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where  £  is  the  location  of  the  source  along  the  x-axis 

G3D  is  a  pulsating  source  translating  at  velocity  U  in  the  x-direction 
dGiD/dy  is  a  pulsating  dipole 
Qj  is  the  strength  of  the  source  distribution 
dj  is  the  strength  of  the  dipole  distribution. 

The  potential  4>j  may  be  taken  to  be  the  corresponding  outer  solution  for  the  inner  radiation  potentials 
4>j,  j  =  1,  ...  ,  6  as  well  as  the  inner  diffraction  potentials  i bjt  j  =  2,  3.  From  symmetry  considera¬ 
tions,  Qj  =  0  for  the  lateral  modes  ( j  =  2,  4,  6)  and  dj  =  0  for  the  vertical  modes  (j  =  1,  3,  5). 

5.1  Relation  to  Near  Field  Singularities 


The  most  direct  approach,  which  will  be  described  first,  is  to  relate  the  strengths  of  this  far  field 
line  distribution  to  the  strengths  of  the  near  field  surface  distribution  calculated  by  using  strip  theory. 
Since  the  cross  section  shape  is  lost  in  the  far  field,  the  singularity  distribution  over  each  cross  section 
may  be  condensed^  to  a  single  resultant  singularity  strength  as  follows.  For  the  symmetric  vertical 
modes  (j  =  l,  3,  5),  tor  which  Eq.  (35a)  applies,  the  resultant  source  strengths  at  the  cross  section 
x  =  o-/£),  is  given  by  the  following  summation 

°V  =  <rj  +  U&j  -  2  £  QjASj,  j  =  1,  3,  5  (48) 

/-I 

where  the  first  equality  is  suggested  by  Eq.  (31) 

iVis  the  number  of  line  segments  modeling  the  +y  half  of  the  cross  section 
As,-  is  the  length  of  the  /th  line  segment. 

For  the  antisymmetric  lateral  modes  (j  =  2.  4,  6),  for  which  Eq.  (35b)  holds,  the  resultant  source 
strength  cr,  is  zero.  However,  the  resultant  dipole  strength  M  is  nonzero 

N  =  +  Ufij  =  2  £  Z)  A s„  j  =  2,  4,  6.  (49) 

/-I  1  -vi 

Newman,  and  Ogilvie  and  Tuck  use  the  method  of  matched  asymptotic  expansions  to  derive  expres¬ 
sions  for  the  outer  singularity  strengths  Qj  and  dj  in  terms  of  the  inner  strengths  cr,-  and  p,  .  The 
essence  of  the  method  is  that  the  inner  and  outer  solutions  are  matched  in  an  overlap  region  r  whose 

distance  from  the  ship  hull  is  large  compared  to  the  beam  B  or  draft  T  and  small  compared  to  the 
length  L 

(B,  T)  «  r  «  L.  (50) 


To  second  order  accuracy  in  e,  where  (B/ L,  T/L)  =  0(e),  the  relation  between  the  outer  dipole 
strengths  dj  and  the  inner  strengths  /j.j  for  the  lateral  modes  is  extremely  simple 


dj  ~  =  My  +  Ujlj,  j  =  2,  4,  6.  (51) 

That  is,  in  slender  body  theory,  there  is  negligible  longitudinal  interaction  between  cross  sections  for 
lateral  motions  and  the  strengths  of  the  inner  and  outer  dipoles  coincide.  (Recall,  however,  that  the 
inner  dipoles  are  two-dimensional  while  the  outer  dipoles  are  three-dimensional.)  The  relation  for  the 
vertical  modes  is  somewhat  more  complex.  In  general,  the  outer  source  strengths  q:  are  related  to  the 
inner  strengths  c r7  by  the  following  integral  equation 


<7/ (*)  ~ 


■-TT  / 


CT; 


~  1 


fL  Qj(0f(x  —  Odi;  =  i Tj(x ),  j  =  1,  3,  5 


(52) 


where  c r j  is  the  complex  conjugate  of  cr 
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fix  -  £;)  is  an  interaction  function  defined  by  Eq.  (4.13)  in  Newman.2  However,  for  the  two  lim¬ 
iting  regions  where 


X  =  0(1) 
X  =  0(e) 


(53a) 

(53b) 


Eq.  (52)  may  be  approximated  by 


Qj  ~  o-j  =  0-9  +  U&j,  j  =  1,  3,  5. 


(54) 


For  these  two  wavelength  regions,  the  inner  and  outer  singularity  strengths  again  coincide  as  in  the  case 
of  the  lateral  modes. 

In  the  present  study,  the  initial  approach  will  use  Eqs.  (51)  and  (54)  to  obtain  the  outer  singular¬ 
ity  strengths.  If  additional  accuracy  is  required,  the  integral  Eq.  (52)  will  be  used  for  the  vertical 
modes. 

5.2  Relation  to  Added  Mass  and  Damping  Forces 

The  second  approach  is  somewhat  more  indirect  than  the  first  approach.  In  this  approach,  the 
computed  forces  acting  on  a  cross  section  are  first  converted  to  a  resultant  inner  singularity  strength 
and  then  Eqs.  (51)  and  (54)  are  used  to  obtain  the  outer  singularity  strengths.  However,  this  approach 
does  have  the  advantage  that  it  would  make  use  of  the  many  previous  calculations  of  the  forces  on  vari¬ 
ous  cross  section  shapes;  see,  for  example,  the  extensive  calculations  by  Porter20  and  the  various  results 
reported  by  Wehausen.21 

One  motivation  for  this  approach  is  that  in  infinite  fluid,  there  exist  the  following  simple  relations 
between  the  dipole  strengths  dy  and  the  added  mass  coefficients  A,)2, 


2D:  (8„  +  Ay/p),  i,  j  -  1,2 

3D:  du  =  ~  (8,j  +  Au/p),  i,  j  =  l,  2,  3 


(55a) 


(55b) 


where  S  is  the  cross  sectional  area 
Fis  the  body  volume 

5 ij  is  the  Kronecker  delta  function  =  0  if  /  ?£■  y'and  =  1  if  /  =  j 
p  is  the  fluid  density. 

In  the  presence  of  the  free  surface,  it  is  the  damping  coefficient  B0  (and  not  the  added  mass 
coefficient  Au)  which  gives  rise  to  far  field  waves,  where  AtJ  and  B,  are  given  by 


(56) 


By  using  Green’s  theorem  on  the  radiation  potential  <f>j  and  its  complex  conjugate  Newman2  obtains 
the  following  relation  between  the  damping  coefficient  and  the  two-dimensional  far  field  wave  ampli¬ 
tude  Aj 


Bjj  =  (pg2/ co2)  \Aj\2,  j  =  1,  2,  3,  4. 


(57) 


By  using  far  field  expressions  for  the  source  and  dipole,  he  further  obtains  the  following  relations 
between  Aj  and  the  resultant  inner  singularity  strengths 


(58) 
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where  a-j  and  fxj  are  respectively  the  resultant  inner  source  and  dipole  strengths  which  give  rise  to 
two-dimensional  waves.  By  combining  Eqs.  (57)  and  (58),  the  desired  relation  between  the  singularity 
strengths  and  the  damping  coefficient  is  as  follows 


a-j  =  2  y/aBjj/p 

My  =  J  JuBjj/p. 


(59a) 

(59b) 


The  above  formulation  applies  for  j  =  l,  2,  3,  4.  Equations  (33)  and  (34)  show  that  <t>5  and  06  are 
entirely  in  terms  of  <a3  and  <*,,  respectively.  In  the  case  of  the  diffraction  potential  <*7,  the  "damping 
coefficient"  may  be  related  to  the  component  of  the  diffraction  force  Ff  which  is  90  degrees  out  of 
phase  with  the  incident  wave.  Ff  is  calculated  in  SMP  by  using  the  Haskind  relation,  Eq.  (38). 

5.3  Kochin  Function  Approach 

The  third  approach  is  to  characterize  the  singularity  distribution  over  the  entire  hull  by  the  Kochin 
function  Hj(k,  9)  which  may  be  represented  in  the  following  two  ways 

(60) 
(61) 

where  9  is  the  direction  of  wave  propagation  with  the  x-axis 
y j  is  the  strength  of  the  singularity  at  rj,  £ 

j  _  ^[kz  +  ik  (£  cos0  +  T}  sin 0)3  (62) 

is  a  wave  elevation  function  at  depth  z  below  the  free  surface  for  a  wave  propagating  in  the  9-direction. 
The  Kochin  function  may  be  used  in  a  variety  of  ways,  including  the  evaluation  of  the  forces  acting  on 
an  oscillating  body  averaged  over  a  period.  Of  principal  interest  in  the  present  study  is  the  fact  that  in 
the  form  given  by  Eqs.  (61)  and  (62)  the  Kochin  function  represents  an  average  of  the  singularity 
strengths  weighted  by  the  wave  elevation  function  f.  In  the  following  section,  an  extremely  simple 
relationship  is  given  for  the  far  field  wave  elevation  in  terms  of  the  Kochin  function. 
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6.  FAR  FIELD  WAVE  ELEVATION 


6.1  Asymptotic  Evaluation  for  Line  of  Singularities 


In  the  first  two  approaches  given  above,  the  far  field  wave  elevation  reduces  to  the  asymptotic 
evalution  of  the  following  equation 


£,e'°"  =  — —  ±  U  .e'"')  =  -1 
g  dt  1  g 


ico  -  U- 


bx 


<t>jeiQit,  j=  1,  ,  7 


(63) 


where  </>7  is  composed  of  the  antisymmetric  potential  and  the  symmetric  potential  i//3.  Since  the 
outer  potentials  are  approximated  by  a  line  of  oscillating  sources  and  dipoles,  Eq.  (47),  the  problem 
essentially  reduces  to  the  asymptotic  evaluation  of  Eq.  (63)  for  a  single  oscillating  source  or  dipole  and 
then  summing  over  all  the  singularities 
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where  M  is  the  number  of  cross  sections  or  strips  and  A Lm  is  the  axial  distance  between  the  mth  and 
(m  —  l)st  strips. 


In  addition  to  the  general  case  U  >  0  shown  in  Eqs.  (63)  and  (64),  the  special  case  U  =  0  is  also 
discussed  at  some  length  since  it  is  considerably  simpler.  One  simplification  for  U  =  0  is  that 
d/dt  —  ia >,  resulting  in  the  following  simpler  expression  for  £j 

*  /< 0  , 

0 - —  <t>J-  (65) 

That  is,  the  value  of  £j  is  simply  an  algebraic  factor  times  the  value  of  (/>,.  The  principal  reason  for  the 
simplicity  of  the  U  =  0  case  is  the  considerable  reduction  in  the  complexity  of  the  expressions  for  the 
source  potentials  as  compared  to  the  U  >  0  case.  Wehausen  and  Laitone22  give  the  following  formulas 
for  the  potentials  G^D  and  G&  for  oscillating  stationary  and  translating  three-dimensional  sources 
placed  at  x  =  a,  y  =  b,  z  =  c 


G!d  =-  +  PV  C  ek’(z~c)  Jo (k'R)  dk' 

r  J°  k  -  k 

^30  —  C  do  C  dk  Fie,  k) 

r  V  \  77  J0  JQ 

Fie,  k)  -  ^'t(z-c)  +  /^~a)cosfli  cos  [k'iy  -  b )  sin  a] 
gk  —  (ca  +  k' U  cos  O)2 

where  only  the  real  parts  of  the  source  potentials  have  been  written 
r  =  VCx  —  a)2  +  {y  —  ~b)2  +  (z  —  c)2 
1/ris  a  Rankine  or  infinite  fluid  source 
PV denotes  a  principal  value  evaluation  of  the  integral 
R  ~  V ix  —  a)2  +  (y  —  6)2 
Jq  is  the  Bessel  function  of  order  0 
rx  =  V (x  -  a)2  +  (y  —  Z>)2  +  (z  +  c)2. 


(66) 


(67) 


The  expression  for  the  integral  term  in  Eq.  (67)  is  a  considerable  abbreviation  of  the  complex  expres¬ 
sion  given  in  [22],  The  actual  integration  interval  is  divided  into  three  parts  to  account  for  the  singular¬ 
ities  arising  due  to  the  various  zeroes  of  the  denominator  of  the  wave  function  F. 


Asymptotic  expressions  for  the  behavior  of  for  large  values  of  R  have  been  obtained  by  a 
number  of  authors  (for  example,  [22],  [23],  [24]) 
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(68) 


where  —  denotes  an  asymptotic  approximation.  The  wave  elevation  £  is  then  obtained  by  using  Eq. 
(65)  and  taking  the  real  part 


£  =  Re 


ho 
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G° 


3D 


—  2tt 
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(69) 


Equation  (69)  is  simply  an  expression  for  ring  waves  with  center  at  the  pulsating  source  and  decaying  at 
the  rate  -Jl/R .  It  should  be  noted  that  even  though  the  pattern  is  circular  around  a  given  source,  the 
pattern  around  an  entire  ship,  which  has  sources  placed  along  the  entire  length  L,  will  not  be  circular. 
Far  field  wave  elevation  for  the  dipole  dG$D/Sy  is  presumably  obtained  by  simply  taking  9/9  v  of  Eq 
(69).* 


*01ver‘  points  out  that  derivatives  of  asymptotic  relations  are  not  always  permissible. 
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Newman  gives  several  references  which  investigate  the  far  field  behavior  of  Gjd  for  large  values 
of  R.  These  include  studies  by  Brard,26  Becker,27  Eggers,28  and  Hanaoka.29  However,  these  references 
do  not  derive  complete  asymptotic  relations.  Instead,  they  largely  focus  on  deriving  certain  general 
features  of  the  far  field  wave  pattern.  All  of  these  studies  as  well  as  the  work  by  Wehausen  and  Lai- 
tone,  show  that  the  behavior  of  the  wave  pattern  is  quite  different,  depending  on  whether  t  is  greater 
or  less  than  1/4  where  r  is  given  by 


Thus,  there  is  the  added  complexity  of  more  than  one  type  of  asymptotic  behavior. 


Two  approaches  will  be  concurrently  pursued  in  the  investigation  of  the  far  field  wave  elevation 
lor  the  general  case  L  >  0.  The  more  straightforward  approach  is  to  continue  a  more  careful  study  of 
[26]  to  [29]  and  additional  references  which  investigate  the  far  field  wave  elevation  due  to  G%.  The 
second  overall  approach  would  be  to  make  use  of  the  extreme  simplicity  of  the  U  =  0  case  given  by 
Eqs.  (66),  (68),  and  (69)  or  reduce  the  extremely  complex  U  >  0  case  given  by  Eq.  (67).  An  obvious 
first  step  is  to  use  Eqs.  (63)  and  (64),  which  contain  the  effect  of  U  in  the  time  derivative,  on  the 
asymptotic  form  for  G$D.  It  may  be  argued  that  this  is  an  inconsistent  approach  in  that  U  is  used  in 
one  part  of  the  calculation  (the  evaluation  of  £)  and  neglected  in  another  part  (the  free  surface  condi¬ 
tion).  This  inconsistent  approach  has  been  used  by  Beck30  and  Lee,  O’Dea,  and  Myers31  (who  use  the 
two-dimensional  source  G%)  to  calculate  the  near  field  wave  pattern  next  to  the  ship  hull.  Also,  SMP 
uses  the  inconsistent  approach  of  neglecting  U  in  the  free  surface  condition  (Eq.  (30))  to  calculate  the 
radiation  potentials  but  then  including  U  in  calculating  the  resultant  forces  (Eq.  (36)),  which  are  rea¬ 
sonably  accurate.  The  accuracy  of  this  inconsistent  approach  for  the  far  field  wave  elevation  will  be 
carefully  investigated. 

A  second,  more  complex,  step  would  be  to  reduce  the  complexity  of  G^D  by  deriving  solutions 
•or  small  values  of  i  in  the  free  surface  condition,  Eq.  (7).  One  approach  would  be  to  neglect  the  term 
involving  U  in  Eq.  (7).  Another  approach  would  be  to  approximate  the  terms  involving  6rand  U2  in 
Eq.  (7)  by  known  functions,  such  as  G®D. 


6.2  Asymptotic  Evaluation  for  Kochin  Function 


Finally,  the  calculation  of  the  far  field  wave  elevation  using  the  Kochin  function  H  (Eqs.  (60) 
and  161))  is  briefly  discussed.  Wehausen  and  Laitone'^  give  the  following  asymptotic  evaluation  of  £• 
derived  by  Kochin,  which  is  restricted  to  the  U  —  0  case 


MR,  e) 


Vi 


K—  Hj(k,  0)  sin 


Rk-f 


.  7=  I,  ....  7 


(71) 


where  Hj  is  the  complex  conjugate  of  H}.  Equation  (71)  shows  that  Cj  is  very  similar  in  form  to  the 
wave  elevation  £  for  a  source  of  unit  strength,  given  in  Eq.  (69).  If  r  =  c  =  0  in  Eq.  (69),  which 
means  that  £  is  evaluated  at  the  free  surface  for  a  source  placed  at  the  free  surface,  the  following 
expression  for  the  ratio  £;/£  results 


f- = ^  vT^(Ar' o)-  (?2) 

The  right  hand  side  of  the  above  equation  may  be  thought  of  as  the  "strength"  of  a  single  resultant 
source-like"  singularity  which  is  equivalent  to  the  weighted  average  of  all  the  sources  on  the  ship  hull. 
Equation  (71)  then  gives  the  wave  elevation  due  to  this  weighted  singularity.  In  the  first  two 
approaches  described  above,  the  contribution  of  each  source  to  the  far  field  wave  is  first  evaluated  and 
the  resultant  wave  elevation  is  then  obtained  as  the  sum  of  all  the  individual  contributions.  The  discus¬ 
sion  following  Eq.  (69),  whjch  states  that  the  resultant  wave  pattern  for  the  entire  ship  is  not  circular,  is 
confirmed  by  the  fact  that  Hj  is  a  function  of  the  direction  6. 
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7.  SUMMARY 


A  somewhat  detailed  description  is  given  of  the  linearized  formulation  for  the  radiation  and 
diffraction  problems.  It  is  pointed  out  that  even  the  linearized  formulation  is  quite  complex.  Various 
simplified  approaches  including  thin  ship,  flat  ship,  and  long  wavelength  theories  are  shown  to  be  either 
unsuccessful  or  inapplicable  to  the  present  interest  of  calculating  the  far  field  wave  pattern. 

The  high  frequency  (low  wavelength)  slender  body  strip  theory,  as  implemented  in  the  widely 
used  DTNSRDC  Ship  Motion  Program  (SMP),  is  identified  as  a  simple  yet  relatively  accurate  approach 
for  calculating  the  flow  field  near  the  ship  hull.  This  theory  neglects  the  forward  speed  effect  in  the 
free  surface  condition  as  well  as  interaction  in  the  longitudinal  direction.  A  unified  strip  theory  which 
accounts  for  longitudinal  interaction  is  mentioned  as  a  possible  later  approach.  Strip  theory  uses  two- 
dimensional  sources  to  solve  the  radiation  problem  due  to  ship  oscillations  in  the  vertical  plane  (surge, 
heave,  and  pitch)  while  dipoles  are  used  for  oscillations  in  the  lateral  plane  (sway,  roll,  and  yaw). 

It  is  pointed  out  that  the  diffraction  problem  generally  leads  to  the  solution  of  a  Helmholtz  equa¬ 
tion  instead  of  Laplace’s  equation  for  the  radiation  problems.  However,  for  wavelengths  which  are 
intermediate  between  ship  beam  (or  draft)  and  ship  length,  Laplace’s  equation  may  again  be  used.  In 
this  case,  the  solution  for  the  diffraction  problem  reduces  to  the  solution  of  a  pair  of  potentials  which 
are  similar  to  the  radiation  potentials  for  sway  and  heave. 

It  is  shown  that  the  near  field  two-dimensional  singularities  are  not  appropriate  for  calculating  the 
tar  field  wave  elevation.  Two  somewhat  similar  approaches  are  indicated  for  converting  the  near  field 
singularity  strengths  to  a  line  of  three-dimensional  far  field  singularity  strengths.  In  one  approach,  it  is 
shown  that  the  strength  of  the  far  field  singularity  at  a  given  longitudinal  station  may  be  approximated 
by  the  resultant  inner  singularity  strength  for  the  cross  section  at  the  given  station.  This  approximation 
is  accurate  to  high  order  for  the  lateral  modes  and  is  reasonably  correct  for  the  vertical  modes  for  long 
wavelengths  of  the  order  of  ship  length  and  short  wavelengths  of  the  order  of  ship  beam  or  draft.  A 
more  complex  integral  equation  approach  may  be  used  if  more  accuracy  is  desired  for  the  vertical 
modes.  In  the  second  approach,  it  is  shown  that  the  far  field  singularity  strength  may  be  obtained  from 
the  calculated  (or  measured)  damping  force  acting  on  the  cross  section. 

Once  the  far  field  singularity  strengths  have  been  determined,  the  far  field  wave  elevation  may  be 
obtained  by  first  calculating  the  contribution  due  to  each  singularity  and  then  summing  over  all  the 
singularities.  It  is  shown  that  the  potential  for  the  stationary  pulsating  source  is  much  simpler  than  the 
corresponding  potential  for  the  pulsating  source  moving  at  forward  speed  U.  In  particular,  extremely 
simple  expressions  have  been  derived  for  the  asymptotic  far  field  behavior  of  the  potential  and  wave 
elevation  for  the  stationary  pulsating  source.  Corresponding  expressions  for  the  dipole  may  be  con¬ 
veniently  obtained  by  taking  lateral  derivatives  B/By  of  the  results  for  the  source.  Asymptotic  expres¬ 
sions  do  not  appear  to  have  been  derived  for  the  pulsating  source  with  U  >  0.  Instead,  most  studies  in 
this  area  (which  are  often  written  in  a  foreign  language)  attempt  only  to  identify  major  features  of  the 
complex  wave  pattern.  It  is  pointed  out  that  the  most  fruitful  approaches  in  this  case  probably  lie  in 
simplifying  the  potential  for  U  >  0  for,  say,  small  values  of  U  to  capture  the  essential  effect  of  U  in 
the  free  surface  condition.  A  particularly  simple  first  step  would  be  to  evaluate  wave  elevation  by 
including  U  in  the  time  derivative  operation  on  the  potential  for  the  stationary  source. 

An  alternate  way  of  obtaining  the  far  field  wave  pattern  is  to  first  calculate  the  Kochin  function 
which  is  a  weighted  average  of  the  singularity  strengths  over  the  ship  hull.  For  the  case  of  a  stationary 
ship,  L  =  U.  the  expression  for  the  resultant  far  field  wave  pattern  is  similar  to  the  expression  for  a  sin¬ 
gle  source  with  "strength"  which  is  proportional  to  the  complex  conjugate  of  the  Kochin  function. 
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Fig.  I  —  Definition  of  coordinate  system  and  six  modes  of  ship  motion 
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Fig.  2  —  Typical  station  offset  distribution 
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